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At redshifts z > 30 neutral hydrogen gas absorbs CMB radiation at the 21cm spin-flip frequency. 
In principle this is observable and a high-precision probe of cosmology. We calculate the linear- 
theory angular power spectrum of this signal and cross-correlation between redshifts on scales much 
larger than the line width. In addition to the well known redshift-distortion and density perturbation 
sources a full linear analysis gives additional contributions to the power spectrum. On small scales 
there is a percent-level linear effect due to perturbations in the 21cm optical depth, and perturbed 
recombination modifies the gas temperature perturbation evolution (and hence spin temperature 
and 21cm power spectrum). On large scales there are several post-Newtonian and velocity effects; 
although negligible on small scales, these additional terms can be significant at I < 100 and can be 
non-zero even when there is no background signal. We also discuss the linear effect of reionization 
re-scattering, which damps the entire spectrum and gives a very small polarization signal on large 
scales. On small scales we also model the significant non-linear effects of evolution and gravitational 
lensing. We include full results for numerical calculation and also various approximate analytic 
results for the power spectrum and evolution of small scale perturbations. 
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The cosmic microwave background (CMB) anisotropies have proved to be a valuable source of information about 
the initial conditions and evolution of the universe. Most current observations measure the CMB temperature and 
polarization assuming an exactly blackbody spectrum. However by looking at the anisotropies as a function of 
frequency vastly more information can be obtained. In addition to the signal from secondary scattering in clusters, 
in principle there is also line-absorption from sources along the line of sight. One of the most interesting of these 
is line absorption due to the 21cm spin-flip transition in neutral hydrogen, giving a low frequency probe of the gas 
• distribution at redshifts 300 > z > 30 pj, 13; 9 • This is sensitive to perturbations on all scales down to the Baryon 
Oh| Jeans' scale, which is orders of magnitude smaller than the photon-damping scale that limits what can be learnt 
directly from the CMB temperature. The 21cm absorption signal therefore potentially contains a huge amount of 
information about small-scale cosmological perturbations. Since absorption signal from redshift z is observed at 
wavelength A = (1 + z)21.106cm [(1 + z)v = 1420.4MHz], the signal can also be studied as a function of observed 
frequency to give tomographic information about the perturbations 0, 0|- Unfortunately observations at many- 
^ \ metre wavelengths are very challenging (see e.g. Refs. 0, Q), but make a useful target for next-but-one generation 
experiments. 

The origin of the dark-age absorption signal is as follows. After recombination there is still a small fraction of free 
electrons. Compton scattering transfers energy between CMB photons and the electrons (and hence the gas), and 
hence keeps the gas temperature close to the CMB temperature until a redshift of z ~ 300. At lower redshifts the 
coupling becomes ineffective and the gas starts to cool adiabatically. Atomic collisions in the gas drive the atomic 
energy levels of the gas towards equilibrium with the gas temperature. The spin temperature defines the relative 
abundance of triplet and singlet hydrogen states, and is driven by collisions towards the gas temperature. Since the 
gas cools faster than the CMB, the spin temperature is below the CMB temperature, and 21cm CMB photons will 
have net absorption by the gas. At redshifts z < 300 an absorption signal may therefore be observable. At redshifts 
z < 30 atomic collisions become very rare, and the spin temperature is driven back towards the CMB temperature 
by interaction with the numerous CMB photons. The absorption signal from the dark ages is therefore limited to 
30 < z 300. At lower redshifts sources of Lyman-a photons and non- linear effects become important, and again 
the spin temperature can depart from the CMB temperature, giving a signal in absorption or emission. 

In this paper we focus on the absorption signal from z > 30 where the physics is well understood and much 
cleaner than the large uncertainties currently surrounding modelling at lower redshifts. We calculate the linear theory 
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angular power spectrum of the 21cm absorption as a function of redshift, including super-horizon scales where post- 
Newtonian effects may be important. We focus on the angular power spectrum Ci(z,z') as this is what is directly 
observable. Many aspects of the physics may be much clearer with a reconstruction of the 3-D power spectrum [f|, 
but converting the observations into such a spectrum is in general non-trivial especially on large scales, and also 
dependent on assumptions about the cosmology. Since the perturbations should be nearly linear at high redshifts 
the statistics should be close to Gaussian, and the angular power spectra should encapsulate most of the statistical 
information in the observation. Our work extends that of Refs. 0, 0, H HI IE EE] by including linear terms due to 
gravitational redshifting, all velocity affects, ionization fraction perturbations, self-absorption, and reionization re- 
scattering. Corrections due to these extra terms are generally quite small, though percent-level effects will be very 
important if high-redshift 21cm is ever going to fulfil its potential for constraining cosmology. We also estimate the 
effect of non-linear evolution, which can be important on small scales even at high redshift, and calculate the effect 
due to gravitational lensing. Although we do not directly consider z < 30 here, many of our results could easily be 
adapted to lower redshifts given a model of the Lyman-a sources and non-linear clustering. 

The CMB temperature anisotropy, sourced at z ~ 10 3 , sees super-horizon perturbations at I < 100. For a 21cm 
signal at z ~ 50 the horizon scale is about ten times larger, corresponding to an angular scale I < 10. One might 
therefore expect post-Newtonian effects to dominate at I < 10. However the small-scale 21cm anisotropy is sourced 
by hydrogen density (and spin temperature and redshift distortion) fluctuations, which grow rapidly towards smaller 
scales. The large scale signal is therefore dominated by fluctuations coming from much smaller scales. Since these 
smaller scales are uncorrelated on large scales, this gives an approximately white-noise 21cm power spectrum on large 
angular scales. This white-noise signal dominates that from super-horizon scales, so the post-Newtonian corrections 
are generally below cosmic variance. In addition there are velocity effects, the most important of which is the dipole in 
the radiation field seen by each hydrogen atom due to its motion with respect to the CMB. This can give non- negligible 
quantitative corrections to the angular power spectrum at I < 100. 

On small scales Thomson scattering of 21cm by the background reionization damps the entire spectrum and also 
induces a very small polarization signal. There are also additional 21cm perturbation sources due to perturbations 
in the 21cm optical depth; for example an overdensity will have a slightly higher optical depth than the background, 
leading to a few-percent suppression in the absorption signal. Additional effects arise indirectly; in particular inclusion 
of ionization fraction perturbations is important for the evolution of gas temperature fluctuations, and can modify 
the spectrum at all angular scales by a couple of percent. 

The approach we adopt is to evolve the Boltzmann equation for the photon distribution function sourced by 
absorption of 21cm radiation by neutral hydrogen. We restrict our attention to a spatially-flat close-to-Friedmann- 
Robertson- Walker cold-dark-matter (CDM) universe. Our results apply equally well with adiabatic or isocurvature 
scalar mode initial conditions, though we only calculate the adiabatic mode spectra explicitly. We note in passing 
that 21cm observations are potentially an excellent way to probe isocurvature modes, especially the compensated 
CDM-baryon mode that cannot be constrained from the CMB temperature [IE [ll[ . We do not consider the tensor 
contribution to the intensity power spectrum as the effect is expected to be well below cosmic variance, though we do 
calculate the tiny tensor-induced polarization signal. We assume no significant particle decay or annihilations, and 
assume no variation of constants or non-Gaussianity, though these can be well constrained from their effect on the 
21cm signal if included pj, El EI El EE E3 • 

We start by deriving the Boltzmann equation for the distribution function in linearized General Relativity. The 
main linear-theory result relating the observable anisotropy on the sky to sources on the absorption surface is given 
in Eq. (fT5|) . On small scales the majority of the terms are negligible, and we give a result accurate on small scales in 
Eq. (|22p . This contains the usual density and spin-temperature fluctuation and redshift-distortion sources, but with 
additional few-percent terms due to the non-zero 21cm optical depth that are often neglected. In Section [HII we then 
derive results for the angular power spectrum in terms of linear-theory transfer functions. To actually calculate the 
power spectrum we need to calculate the sources, so in Section ITVl we give results for calculating the background and 
perturbed densities and temperatures. In Section|V]we describe the qualitative shape of the power spectrum, and give 
approximate semi-analytic results for the large and small scale power spectrum. We then quantify the importance of 
the various effects in Section [VII where we calculate the 21cm intensity power spectrum numerically The very small 
large-scale polarization signal is calculated in Section IVlII Since the 21cm power spectrum probes small scales, non- 
linear evolution can in fact be important at the many-percent level even at redshift z ~ 50. We give an approximate 
estimate of the effect in Section IVlIIl and also give an accurate result for the lensed power spectrum in Section [DTI A 
full non-linear analysis is beyond the scope of this paper, but we briefly discuss other sources of non-linear power in 
Section [Xj We finally conclude in Section IXll In a series of appendices we give approximate results for the evolution 
of small-scale baryon perturbations, general-gauge results for numerical calculation, equations for the evolution of the 
ionization fraction perturbations, useful results for integrating spherical Bessel functions, and third-order perturbation 
theory results for the non-linear CDM density and velocity power spectra. 

The Boltzmann approach to calculating the line emission angular power spectrum that we consider here is related 
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to the angular power spectrum of source number counts; we discuss this in a companion paper [18j |. 

II. BOLTZMANN EQUATION 

At high redshift, perturbations should be close to linear, and we assume there are no sources of Lyman-a photons. 
During the dark ages the atomic collision time is comparable to the CMB photon interaction time, and for accurate 
results the full distribution of spin and velocity states must be accounted for [l9j . To simplify our analysis we 
neglect this complication, focussing on new effects that arise from a full linear perturbation analysis when the spin 
temperature is independent of atomic velocity. The spin temperature during the dark ages in then governed by 
2 lcm- interaction with CMB photons and atomic collisions. Since the collision rates are not known very accurately, 
and the ionization fraction after recombination is somewhat uncertain, the precision of our calculation is currently 
limited anyway. We also assume the background CMB temperature is exactly blackbody, neglecting any effects due 
to non-21cm distortions, and assume that atomic angular momenta are isotropically distributed. 

We employ linearized standard General Relativity, working in the conformal Newtonian gauge with metric 

ds 2 = a 2 {r])[(l + 2i>)dtf - (1 - 2(j))5i j dx i dx j }. (1) 

Except when quoting a few results relevant for calculation of numerical answers we use natural units with c = 1. We 
take a velocity field u a to be along d v so that u M = a _1 (l — i/j)Sq an< ^ u n ~ a (l + VO^o • This velocity field is the 
zeroth element of an orthonormal tetrad which we take to be (Ao) a = u a and Xi = a _1 (l + cj))di. Decomposing a 
photon wavevector k a = dx a /d\ into a direction e a and frequency k ■ u = e/a, we have 

|p(l + + ^)e, g=a- a e(l-V0, (2) 

where the three- vector e comprises the spatial components of the propagation direction on the spatial triad Xi. 

The number density of neutral hydrogen atoms is nni = no + ni, where the density in the ground state (degeneracy 
1) is no, and the density in the upper triplet state (degeneracy 3) is ri\. We assume the spin temperature T s is 
dependent only on time and position, defined by rii/no = 3e~ T */ Ts where T+ = h v v 2 \jkB ~ 0.068 K and v^i is the 
constant 21cm frequency. In the rest frame of the gas the net number of 21cm photons emitted per unit volume in 
proper time <h~ g within energy dE within solid angle dfi is 

dn 21 = — [(ni - 3n )AC + n x ] A 10 <S>(E - E 21 )dT g dEdfl, (3) 

47T 

where E 21 corresponds to the 21cm frequency and J\f v is the photon phase space density controlling stimulated 
emission. The line profile Q(E — E 2 i) is defined so that J dEQ(E — E 2 \) = 1. The spontaneous emission rate 
A\q = 2Trav 21 hp/ (3c 4 ml) « 2.869 x 10 _15 s _1 [2(J, corresponding to a spontaneous decay time of ~ 10 7 years and 
CMB photon interaction time ~ T*/(T 7 Aiq) (about 10 4 years at z = 30). The form of the equation follows from 
considering detailed balance in equilibrium, in which there is no net production of photons at any temperature. 

We shall only make an accurate calculation on scales larger than the line width. In this approximation we model 
the source emission as monochromatic, so that $(E — E 2 i) = 8{E — E 2 i). The thermal line width corresponds to 
scales with wavenumber of a few hundred Mpc -1 0, [HI- This unavoidability suppresses observed power on very 
small scales regardless of the observational bandwidth. On these scales power is also suppressed due to significant 
baryon pressure at earlier times, as discussed further below. Observational bandwidth can be accounted for simply 
by integrating our final result over a frequency window function. 

We model the radiation field as a CMB blackbody J\fp = (exp(h p v / ksT^) — l)^ 1 plus a term due to 21cm emission 
J\ff. At the temperatures of interest where -C T 7 a good approximation is Afp(v = v 2 \) — T~ i jT i , so that 
M v = Ty/T* +Nf{v = v 2 \). Usually the photon temperature T 7 is taken to be isotropic, but here we are interested in 
corrections and so allow for its angular variation (e.g. due to the dipole in the rest frame of the atom). The gas-frame 
temperature is given by 

TW(e ) x,7 7 )=T 7 ('7)[l + 0(e,x,77)-e-v(x,7 ? )], (4) 

where O is the temperature perturbation and v the gas (baryon) velocity relative to u a on the Xi triad. Anisotropy 
in the radiation field may result in an anisotropic distribution of the atomic triplet states, which would significantly 
complicate our analysis. This will not be an issue if atomic collisions isotropize the distribution rapidly compared to 
the photon interaction time. However during the dark ages the collision and photon interaction times are actually 
similar, so this may not be a safe assumption. Nonetheless, because the interaction is parity invariant, odd multipoles 
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of the radiation field will not affect the triplet distribution; in particular the dominant dipole term leaves an isotropic 
distribution unchanged. Higher order CMB anisotropics will drive an anisotropy in the triplet distribution, however 
their relative amplitude is ~ 1CP 4 so their contribution is very small. Also the scattering time for a photon from 
the anisotropic part of the distribution will be longer than the collision time, so the distribution should in any case 
be randomized effectively. Similar comments apply to the effect of polarized radiation, so our approximation of an 
isotropic triplet distribution should be accurate. 

If the gas 4- velocity is the rest frame energy is given by E^ = k a u a g = e(l — e • v)/a, and an interval dA along 
the photon path corresponds to a proper time dr g = u a g dx a — k a itgd\. The Boltzmann equation for the evolution of 
the distribution function / (number density of photons fdQE 2 dE/c 3 = 2NfdQv 2 dv / c 3 ) due to 21cm interaction is 
then 
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(5) 



where uri is the number density of neutral hydrogen and we used the good approximation T* <C T s . Wc defined 
fractional perturbations in a quantity X as Ax = (6X)/X and denoted background quantities by an over-bar. Since 
the baryon pressure is very low, and the ionization fraction in the dark ages is small, we may take A^i — Ah = Af,, 
though the baryon perturbation A;, can differ significantly from the CDM perturbation A c . 

Although we do not model 21cm emission from reionization in detail here, we do include re-scattering of emission 
from higher redshift by the background electron density as this affects the power spectrum from the dark ages. For 
the moment we neglect polarization and discuss this later. The Thomson scattering contribution to the Boltzmann 
equation is then 
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where F and f% are the monopole and quadrupole parts of /. In the first line here, the tildes denote quantities in the 
gas frame evaluated on the Lorentz-boosted tetrad X^. 

The background equation does not depend on the Thomson scattering term. Defining 
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the background equation becomes 
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where the background optical depth to 21cm is defined by 
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Here e = a(rj e )E2i, 9{x) is the Heaviside function and A denotes the observation point; a subscript e denotes the 
quantity is evaluated at time 7j e [and additionally position + 11(77,4 — rj e ) for perturbed quantities along a line of 



5 



sight n]. The optical depth r e is quite small, typically 1-4% over the epoch of most interest (see Fig. [2] below). The 
time derivative is given by 



f = a,T s 5(e/a - E 2 i) = T e 5(rj - rj e ), 



(10) 



which defines the optical depth source t s (jj) in analogy with p s 
The background solution is then given by the integral of Eq. 
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where TL is the conformal Hubble parameter and /(e) is the value of f(r],e) at rj > r\ t . To first order in r e one can 
use f(r], e) = f(e)9(r) — r] € ), however the full form given above must be used to get results correct to higher order in 
T e . Predictions for the 21cm power spectrum arc often quoted in terms of the brightness temperature today, given by 
7b = E b s hpf /2&b- The isotropic brightness today due to background emission is therefore 



(12) 



During the dark ages the spin temperature is below the CMB temperature, so TJ, is negative corresponding to net 
absorption. 

To calculate the perturbation to the distribution function we define the monopole source 
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The total perturbed Boltzmann equation then becomes 
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Here + denotes the gauge-invariant temperature anisotropy sources with I > 2, and v 7 is the velocity (i.e. dipole) 
of the photon distribution. 

To solve the perturbed Boltzmann equation we use the time component of the geodesic equation, which reduces to 
an equation for the evolution of the comoving energy along the line of sight: 
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(15) 



where overdots denote conformal-time partial derivatives. We parameterize the distribution function in the Newtonian 
gauge as /(?7,x,e,e), in which case the Boltzmann equation becomes 
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where r c is the Thomson scattering optical depth. Noting that ed e = 7i e 1 d r)e , we have from Eq. (jTTJ) 
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which defines an additional time- independent quantity /,i n e(e). Substituting into the Boltzmann equation and into- 
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grate formally along the background line of sight gives the final result 
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where n = — e|^ and e < au-E^i (i.e. observed energy strictly less than E^i), and <5/o is the monopole perturbation. We 
have assumed that the CMB anisotropics can be well observed at higher frequencies, and the first-order perturbation 
from last scattering subtracted off the 21cm map. The 21cm brightness fluctuation calculated here is then that of the 
difference map without the blackbody contribution. We have also assumed there is no overlap between reionization 
and the 21cm absorption. 

The first term in Eq. (fT5|) is the usual monopole source, with additional terms if) + n • v reflecting the additional 
emission due to the difference between proper time in the gas frame and the interval d?7 along the line of sight. Then 
there is the effect from the CMB dipole seen in the gas frame, v 7 — v, plus higher multipole contributions to the 
temperature seen by the source, 0+. The remaining terms on the first line of Eq. (|18|) describe the local effect of 
gravitational and Doppler redshifting on the relation between an observed redshift interval Az and the Ar/ along the 
line of sight. The main such redshift-distortion effect comes from the radial gradient contribution to n • (dv/drj). 
The term multiplying f t i ne (e) is just HSr] where Sr] is the perturbation to the conformal time for a fixed (gas-frame) 
redshift surface; it has the usual Doppler, Sachs- Wolfe and integrated Sachs- Wolfe contributions familiar from CMB 
studies. The first term on the third line describes the perturbation to the 21cm optical depth: for small r e the term 
is proportional to r e /2, corresponding to 21cm photons on average seeing half of the perturbation along their line of 
sight. 

To understand the way in which reionization enters Eq. (|18[) , consider reionization approximated by a delta- function 
visibility function at time rj rc with optical depth r c . Dropping self-scattering terms (taking r e — > 0) we then have 
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The first two sets of terms are the Sf without reionization mulitplied by the fraction, e~ Tc , of 21cm photons that are 
not rescattered. The third set contains the effective 8r] for those photons that do rescatter (i.e. only the common part 
that is accrued after reionization) weighted by the fraction, 1 — e~ Tc , that scatter. Finally, the fourth terms arise 
from in-scattering at reionization and represents an average of the source functions on the electrons' 21cm surface. 
For perturbation modes with fc(?7 rc — r) e ) 3> 1, the dominant contribution to the 21cm monopole at reionization is 
<5/o(e)|rc ~ /.ine(e)V'|rc since the source terms on the electrons' 21cm surface average to zero. Using this in Eq. (|19| . 
we see that for such modes reionization damps the 21cm anisotropies by e~ Tc . For modes with fc(r7 re — ?7 e ) -C 1, 
reionization has essentially no effect since scattering out of the line of sight is balanced by in-scattering. However, 
the contribution from such modes (which were necessarily outside the horizon at r) e ) to the 21cm anisotropy at any 
multipole I is now small since there is considerably more power in modes with larger k (see Sec. [V}. The net effect 
is that reionization should suppress the 21cm anisotropies by e _Tc on all scales, unlike the CMB where there is no 
suppression at large 

On small scales (well inside the horizon) most of the terms are completely negligible. Defining r T = r c e _Te /(l — e~ Tc ) 
Eq. (fl8|) is approximated by 
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(20) 
(21) 
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where x = Va — V is the conformal distance along the line of sight. The result in the first line should be very accurate on 
small scales. In the second line we made the approximation that r e ~ 0, r T ~ 1 to recover the standard approximation 
that is only accurate to C(r £ ) (percent-level). Neglecting small photon perturbations, Eq. ([20|) can also be written as 
an expression for the brightness perturbation today 
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To consider the impact of the additional terms on large scales we next derive an expression for the power spectrum 
for numerical calculation. 



III. ANGULAR POWER SPECTRUM 



For numerical work one can expand into multipoles and harmonics. We use 
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for the ith species, and similarly for the temperature multipoles, giving 
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for I > 1, where a prime denotes a derivative with respect to the argument. The 0; are the angular moments of the 
Fourier expansion of the CMB temperature anisotropy and are defined analogously to Fi. The last term is small, of the 
order of the I > 2 CMB temperature fluctuation. Note that v + Hv — kip is zero in the absence of Thomson scattering 
or baryon pressure effects. Equation (|25p can be integrated over a given frequency window function (determined by 
the observation) to determine the actual observed power. If desired one can integrate by parts so that the result 
depends only on ji and derivatives of the window function and sources. We perform our numerical calculations this 
way in the synchronous gauge: the equations in a general gauge are given in Appendix iBl 
The angular power spectrum is given by 



C ; (z,z')-4tt / dlnkV x (k)F l {k,z)F l (k,z'), 



(26) 



where V x is the power spectrum of the primordial curvature perturbation x an d Fi(k, z) is the distribution function 
multipole transfer function to redshift z — a^I^i/e — 1> 1 - e - z ) — Fi(t]a, E21/O- + z ), k) for unit initial curvature 

perturbation. 

To calculate the sources for the line-of-sight integral we need to compute the perturbations in the spin and gas 
temperatures, and, for reionization, the evolution of the low multipoles of the distribution function. We consider these 
next. 
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FIG. 1: Evolution of the interaction times for H-H, H-e, H-p and H-photon spin-coupling processes, and how this influences 
the spin temperature T s relative to the background CMB and gas temperatures. At high temperatures the H-H collision time 
is short and collisions couple T„ to the gas temperature T g ; at lower redshifts the gas is diffuse and CMB photon interactions 
drive T 3 to the CMB temperature T 7 . This figure assumes purely linear evolution and no Lyman-alpha coupling; in reality 
non-linear effects are likely to change the result at z < 30. 



IV. EVOLUTION 



The evolution of the spin temperature is determined from the evolution equation for the states of a fixed number of 
atoms Nhi = Nq + N\ in the gas rest frame. If we crudely assume that recombinations are to the singlet and triplet 
state in the ratio 1:3 we have 

^ = -N (C 01 + 3A 10 ACo) + Wi(Ci + A 10 (l+M v0 )) - + — , (27) 

where r is the gas proper time, A/^o the monopole part of M v in the gas rest-frame integrated over the line profile 
(evaluated in Appendix [Cj) . and the ionization fraction is x e = n e /riH = n e /(jiHi + n e ) (we are neglecting molecular 
hydrogen and assume all helium is neutral) . Here the collision term is Cio = k^ q h tihi + n\Qn e + K^np and 
Coi = 3Cioe~ T */ Ts where T g is the gas temperature. The rates k\ for hydrogen-hydrogen, hydrogen-electron and 
hydrogen-proton collisional coupling are taken from Refs. [H, [2l[ (fit using a cubic splines and 5th order polynomial 
in the logs), and not known to an accuracy of better than a few percent. As shown in Fig. [1] the hydrogen- hydrogen 
term dominates because of the small dark-age ionization fraction (x e ~ 2 x 10~ 4 ), and proton- hydrogen rates are 
suppressed relative to electron-hydrogen rates by a factor of ~ (m e /TO p ) 1,/2 due to the lower proton velocity, except 
at low redshifts where the proton cross-section is significantly higher plj . We include the small correction from 
electron- hydrogen collisions, but neglect the proton- hydrogen term. 

At redshift z ~ 70 the total collisional and photon interaction rates are about equal, Cio ~ T^Aiq/T*, with 
Cio ~ 10 _11 s _1 (corresponding to a collisional coupling time ~ 4000 years). At lower redshifts the gas becomes diffuse 
and collisions become less effective at coupling the spin temperature to the gas temperature. Defining (3 S = 1/T S 
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FIG. 2: The background 21cm brightness Tb, optical depth r e , and [1(1 + l)Ci/2n} 1 ^ 2 at I = 10 4 as a function of source redshift. 
The dashed line shows the result for Tb neglecting the second term in Eq. pip due to the effect of absorption on the ambient 
blackbody spectrum. 




FIG. 3: Evolution of the fractional baryon, matter and spin temperature perturbations as a fraction of the CDM density 
perturbation. The left figure is for a k = O.lMpc -1 mode, the right figure shows the effect of baryon pressure at k = 500Mpc _1 . 
The dotted lines show the equivalent results neglecting ionization fraction perturbations. In both cases the baryon perturbation 
is significantly less than the CDM perturbation at all relevant redshifts. There is no reionization. 
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(etc.), the evolution of the spin temperature is determined by 

d(3 s (3 S dx 



dr 1 — x e dr 



= 4 [(J3 g - p,)C 10 + & (1 - p, [T 7 + T b0 ] ) A 



10 J 



(28) 



where T 7 + T&o = T*N v o is the (perturbed) monopole brightness temperature. For scales large compared to the 
narrow line profile, we have 



-t v (T s 



T 7 (l + A TT ) + ^(T s -r 7 ) 



i + a s + i> + — 



to first-order in the small optical depth r n = r £rj (see Appendix |C|) . The quantity 



32nk B T s V a u a g 



(29) 



(30) 



is the perturbed optical depth to 21cm, where V a u g is the volume expansion rate of the gas. Note that T s r n is 
independent of the (perturbed) spin temperature. Since t v is small, and the relevant non-perturbative equation in 
t v cannot be easily solved, we make this first-order approximation below. In the epoch before reionization the spin 
temperature and ionization fraction only vary on Hubble time scales. The coupling time is short compared to the 
Hubble time, so the spin temperature is determined to very good accuracy by equilibrium, with the left hand side of 
Eq. being zero, 



c 10 ?yr 7 + a 

C w n/T g + A 



10 



10 



1. 



-T s T n A w 



C W T+/T g + A w C 10 n/T 7 + A 



10 



(31) 



The spin temperature varies between T 7 and T g depending on whether the radiation or collision terms dominate; see 
Fig. [TJ The second term in Eq. (|3"Tj) due to the finite 21cm optical depth is generally very small, giving a correction to 
the spin temperature of less than half a percent, and to T 7 — T s of at most about one percent. The small effect on the 
brightness is shown in Fig.fS] This is smaller than the correction due to our assumption of a single velocity- independent 
spin temperature [l9j . 

The perturbations to the spin temperature are determined by 



(i? 7 - 
Ami 



R g )5C w + C w {R g A Tg - i? 7 A T7 ) + I Tni 4 10 C 10 ^jr T ' 



-iv-v 



A Cl0 + 2A Tt (C 10 -R 7 - 1) - 2R 1 5C_ 
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where Rr l = (C w + A w Ti/T*) and 
^Cio = A^Cio 



[1 X e )KiQ -f 



dlnT„ 



dlnT„ 



T —T 



n H A Tg + (kJ* 



HO 



An 



A T . 



)x e n H A Xe 



(32) 



The background ionization fraction x e is taken from RECFAST [22| . The effect of the optical depth term in Eq. ([32 
is only about 1% on the angular power spectrum. 

Assuming purely Compton cooling, the background gas temperature evolves according to [13, [23| 



T a + 2HT a 



■iaCFTP-yXe 



3m e c(l + / Hc + x e ) 



{T g -f, 



1) 1 



(33) 



where ot is the Thomson scattering cross-section, /jj e = nHc/^H and we have ignored the very small effect of 21cm 
radiation. The perturbations evolve with 



a t =24>-hv ^ihh — 

3 3m e c(l + /ho + x e ) 



l- f I)<|4A T ,+^ 



1 + Z e /(1 + / He ) 



T 

± 9 



(34) 



where we neglected helium fraction perturbations. Note that although the direct contribution of A x<! in Eq. 
is small, the indirect effect on the evolution of At can be significant. The equations for approximately calculating 
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FIG. 4: Transfer function for monopole source at redshift z — 50 given unit initial curvature perturbation, compared to 
other relevant perturbations. The perturbations are numerically evaluated in the synchronous gauge. The Newtonian-gauge 
fluctuations equal those in the synchronous gauge well inside the horizon (k 10 -3 Mpc^ 1 ). On large scales, the Newtonian- 
gauge A;, flattens out at ~ 6/5 times the primordial curvature perturbation, as shown by the dotted curve. The spin temperature 
perturbation is negative at z = 50 (see Fig. [3}. 



A Xe are given in Appendix |Dj An overdensity has positive At but recombines more fully than the background and 
hence has negative A Xe ; the additional effect of the ionization fraction perturbation is therefore to reduce the coupling 
to the CMB and hence slightly decrease the spin temperature. Typical transfer functions for the perturbations are 
shown in Fig. |4] in the synchronous gauge that we use for numerical work. The Newtonian-gauge functions differ on 
super-horizon scales (fc < 10~ 3 Mpc -1 at z = 50). For example, the CDM transfer function flattens to « 6/5 on large 
scales. 

Assuming an ideal gas, the gas pressure perturbation Sp/p g — c 2 s A g is given by 



2 A _ W 



A, 



-2-9- ([1 + x e ](l - F He ) + F He mp/mHe) (A s + A Tg 



(35) 



where /i is the mean particle mass and Yh c is the mass fraction in helium. This result must be used on scales 
where the baryon pressure is important [9(. Note that on these scales there may also be additional effects due to 
CDM decoupling [24| that we neglect here. The evolution of the gas and spin temperature perturbations is shown in 
Figure[3]for two different scales, along with the relative evolution of the baryon and CDM perturbations. Again, these 
are in the synchronous gauge but they are very similar to the Newtonian-gauge perturbations since both wavelengths 
are sub-horizon for the redshift range plotted. 

To evaluate the sources for reionization Thomson scattering we need to evolve the multipole equation 



7 C 

Ft + ^-pY [(I + l)F l+1 - + 5 lo edj4> + -fedJW 



ap s 



6{e/a-E 21 ) 



f [F t + 6io(A H i - A Ta + ^)f] ~ fa- (d e [ej) + fedj) . (36) 



Neglecting the self- absorption terms involving f this is straightforward to propagate numerically after integration 
over a window function in frequency to remove the delta-functions, as discussed further in Appendix [Bj Other 
perturbed quantities and the photon temperature and polarization multipoles evolve according to standard results as 
implemented in the numerical codes CMB FAST and CAMB [H [H, H3, tH • 
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V. APPROXIMATE RESULTS 



The general form of the equal-redshift angular power spectra can easily be understood. On super-horizon scales the 
potential 4> is close to scale invariant and constant, so (from the Poisson equation in the comoving gauge) A c ~ k 2 /H 2 (j>. 
Hence on entering the horizon (k ~ TL) the matter perturbations are of order of the potential, A c ~ (f>. During radiation 
domination photon pressure prevents gravitational collapse: the perturbations only grow logarithmically. As a result 
the spectrum of A c is approximately scale invariant on sub-horizon scales at matter-radiation equality, though super- 
horizon scales still have A c ~ k 2 /TL 2 <j). During matter domination the potential remains constant and A c grows at 
the same rate independent of k on all scales where baryon pressure can be neglected. The result is that an initially 
scale-invariant potential gives comoving-gauge CDM perturbations with an amplitude scaling as k 2 on large scales 
and growing logarithmically with k on small scales. The Newtonian-gauge CDM perturbation equals that in the 
comoving-gauge well inside the horizon, but on larger scales the Newtonian-gauge A c « —2(f) in matter domination. 

For super-horizon modes, the dominant 21cm sources in the Newtonian gauge are 

Sf( VA , x A , e, n) « /(e) (A s + - /, lne (e)V, (37) 

where we have ignored reionization and the effect of non-zero 21cm optical depth. The source A s scales with the 
hydrogen perturbation Ajjj. After recombination almost all the atoms are neutral and Ahi = A^. On scales above 
the baryon sound horizon at recombination the baryons fall into the CDM potential wells on a Hubble time scale, so 
At, evolves to follow A c . Note that although Af, qualitatively follows A c well after recombination, the difference can 
be tens of percent on all sub-horizon scales at high redshift. On very small scales growth of Ah is suppressed once the 
perturbation reaches pressure support. The Newtonian-gauge 21cm sources thus have a scale-invariant amplitude on 
super-horizon scales, scale as k 2 for sub-horizon modes that entered the horizon after matter-radiation equality, are 
growing logarithmically on small scales, and are suppressed on very small scales where c 2 k 2 /TL 2 > 1 at recombination. 
This behaviour can be seen in Fig. 0J 

On small scales the fractional source fluctuation is of the order of the CDM density perturbation, A s ~ (k/TL) 2 (f), 
and the velocity is given by v ~ (TL/k)A c . For small scales with k ^> TL the line-of-sight result in Eq. (|25[) is therefore 
dominated by the monopole and redshift distortions effects, giving the usual approximation for the 21cm source when 
we neglect self- absorption effects (take r £ =0): 



(38) 



Fi(r) A ,x A ,e,k) w e T <=/(e) 
The fractional angular power spectrum for one redshift shell is then 
CHz, z)«47re- 2T « / d\nk\v As (k, z)[ 3l (k x {z))] 2 + 2V V A,(k,z)ji(k X (zM(k X (z)) + V vv (k, z)[f/ (k x (z))} 2 \ , (39) 



where Va b , V v , and V v a b are the power spectra of A s , kv/TL, and their cross-correlation 1 . For k X several oscillations 
larger than I, and smooth power spectra, we can replace the products of the rapidly oscillation Bcsscl functions with 
their approximate smooth averages from Appendix [El 

One might expect the low-Z 21cm fluctuations to be dominated by the super-horizon scale, post-Newtonian fluctu- 
ations at scale k = l/ X . This is not correct since small-scale fluctuations, with amplitude growing rapidly with scale 
as k 2 , couple to I < k X through the oscillatory tails of the spherical Bessel functions. The net effect is that, for all 
I and a single source plane, the dominant contribution is from modes that are inside the horizon. At low I, the C\ 
integral of Eq. (|26|) is therefore dominated by much smaller scales, with k 3> l/ X - I n this limit the Bessel functions 
can be approximated with the asymptotic result ji(k X ) ~ cos [k X — (/ + l)7r/2] /(k X ). Since the power spectrum is 
quite smooth we can then average over oscillations using (\ji{k X )\ 2 ) ~ l/[2(k X ) 2 }. We can similarly remove oscillating 
terms in the second derivative term. Hence on large scales, keeping only the monopole source and redshift distortions 
and assuming a narrow redshift window, the dimensionless fractional power spectrum is 

G(z,z)^4. e ~ 2 -|dlnfc^|^, (40) 

where V s (k, z) is the power spectrum of A s — kv /TL at redshift z (conformal distance X (z)). Since this is independent 
of I it corresponds to a white noise spectrum. The intrinsic fluctuations on a scale k — l/ X are hidden beneath random 



1 We define power spectra so that (A(x) 2 ) = J dlnfc'P(fc). 
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variations in the large scale distribution of much smaller perturbations. This is a generic feature of the 21cm angular 
power spectrum that is also true after the dark ages [f|. 

Scales inside the horizon at matter-radiation equality have an approximately scale-invariant spectrum (grow loga- 
rithmically with k). This causes the 21cm power spectrum to flatten out. First consider the case where the window 
function is sharp, Ax/x ^ 1/h so the source is from a single redshift. For monopole and velocity sources with 
power-law spectra, the result for Ci from Eq. ()39|) can be obtained analytically using a result for integrating products 
of spherical Bessel functions quoted in Appendix [E] In particular, taking the power spectrum to be approximately 
constant on small scales, we can approximate the dimensionless fractional power spectrum as 



2tt 



-Ci{z,z) — e" 



PA.(7rf/2x(*), z) - \v vA , {3nl/4 X (z),z) + lv v (15nl/16x(*),z) 

6 5 



(41) 



The numerical factors are consistent with the angular average of Eq. (|21|) . Since the small scale spectrum actually 
grows logarithmically, the power spectra in Eq. I|4ip are approximated by their values at the mean position of the 
corresponding window function. The Bessel functions are skewed to k > l/ X so that the mean of a [ji(r)] 2 /r window is 
at r = 1(1 + l)n/(2l + 1) ~ ttI/2. The velocity and cross-power window functions probe even smaller scales ~ 157rZ /16 
and ~ 37rZ/4 respectively. Eq. (|4Tj) is accurate at the 10%-level from the end of the baryon oscillations to the 
baryon damping scale at I ~ 10 6 . The power is overestimated because the underlying power spectra are only growing 
logarithmically rather than linearly. In the (crude) approximation that A s ~ A& ~ A c so that kv/Ti. ~ kv c /Ti. ~ — A s , 
and taking V(k) to be constant we have 

l l±^ Cl (z,z)^ 2 ^e^V As (k,z). (42) 

Though not very accurate this result shows the importance of redshift distortions: it is ~ 28/15 ~ 1.87 times larger 
than the equivalent approximate result neglecting them. 

On scales where the wavelength is much smaller than the redshift bin width, I > x/Ax, redshift distortion effects 
average out and we can instead use the Limber approximation: 

d(z,z) « e~ 2T ^£' xdxW(x) 2 V As (k = l/x(z); Vo ~ *(«))> (43) 

where x* is the far end of the window function W(x)- Since the source power spectrum is nearly scale invariant, 
Ci therefore scales approximately as l/l 3 . In other words l 2 Ci is approximately constant with an additional l/l 
suppression due to line-of-sight averaging through the window. If W is a Gaussian of width c, the line-of-sight 
averaging causes the overall amplitude to scales approximately as 1/er. The Limber approximation can in fact be used 
for numerical calculation on scales with I 3> x/^X where it becomes accurate. 

One very small scales where modes are inside the baryon sound horizon at recombination, kc s (?7*)?7* > 1, the 
baryon pressure becomes important and Ah differs significantly from A c even at late times. We discuss approximate 
analytic solutions for the evolution of A& in Appendix [A] In the small scale limit where k 2 c 2 3> Tt 2 the pressure 
and gravitational forces approximately balance, and Ah/A c ~ 7i 2 /(k 2 c 2 ). Since A c is roughly scale invariant, and 
neglecting the effect of the baryon pressure on the CDM evolution, this implies l 2 Ci cx l/l 4 , giving a characteristic 
sharp fall-off in power on very small scales (there is an additional power of l/l for window- function line-of-sight 
averaging). Note the observations on such small scales are unlikely to be possible at high redshift over most of the sky 
due to scattering by turbulent galactic and solar-system plasma [29( . The non-zero linewidth also becomes important 
one these scales, so our results in the approximation of a monochromatic source will overestimate the power. 



VI. NUMERICAL RESULTS 



For our numerical results we assume a standard concordance flat adiabatic CDM (Vl c h 2 = 0.104) model with a 
constant primordial spectral index (n s = 0.95) and optical depth to Thomson scattering t = 0.09. We take baryon 
density Vl^ti 2 — 0.022, Hubble parameter 73kms _1 Mpc _1 and initial curvature perturbation power on 0.05Mpc _1 
scales A s = 2.04 x 10 -9 . We neglect the neutrino masses, which do not have a large effect for source planes at high 
enough redshift that the neutrinos are still relativistic. We take our window function to be a Gaussian of width Av 
in frequency over the observed brightness. 

Fig. O shows the effect on the angular power spectrum of allowing for self-absorption and ionization fraction 
perturbations. The finite non-zero optical depth lowers the amplitude both because the background signal is lower 
and because the optical depth from an overdensity is higher than from the background; each effect is about r c /2 ~ 1.5% 
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FIG. 5: The effect of perturbations to the optical depth and ionization fraction on the 21cm power spectrum at z = 50 with a 
sharp window function. The solid line shows our main result, the dashed-dotted line is the larger result using the a zeroth-order 
expansion in r 6 , the dashed line is the lower result if the effect of ionization fraction perturbations on the gas temperature 
evolution is neglected. The fractional change in the spectrum is roughly the same on all scales where baryon pressure is 
negligible. 



in amplitude, giving an overall suppression of ~ 6% in power. Ionization fraction perturbations increase the power 
by about 2% on all scales: overdensities recombine more fully and hence have less Compton-coupling to the CMB 
and hence lower spin temperature. These effects are important on all scales and should be included in any accurate 
calculation, though note that they are comparable to others effects that we have neglected because of the simple 
velocity-independent spin-temperature approximation (see Ref. [HI). 

On small scales the post-Newtonian and extra velocity terms can be neglected to good accuracy. On large scales 
they can be more important. In Fig. [5] we show the contribution of the auto-variance of various terms to the total 
large-scale power spectrum at a given redshift. As expected the dominant contributions are still from the monopole 
source fluctuations and redshift distortions 2 . Except on very large scales the next most important term is from the 
CMB dipole in the rest frame of the gas, which gives percent level contribution at Z < 50 for the redshift shown 
here, but is negligible on much smaller scales 3 . At (I < 100) the contribution from the potential and other velocity 
terms are also not entirely negligible. The contributions from the CMB temperature anisotropy above the dipole, 
and reionization re-scattering sources are completely negligible on all scales. At lower redshifts the background signal 
becomes smaller, and the relative importance of the terms changes. The background signal depends on T s — T 7 , 
but some of the perturbation sources depend only on the 21cm optical depth and are non-zero even when the spin 
temperature is equal to the CMB temperature. As an extreme example, Fig. [7] shows the relatively large contribution 
from the photon-baryon dipole at z = 20 on large scales if there were no additional sources from non-linear structures. 

Note that just because something does not show up in a narrow frequency window auto-power spectrum at a given 
redshift does not mean that it is necessarily negligible. The correlation between source planes at a given I falls off very 
rapidly once the plane separation is greater than characteristic perturbation size Extra information may therefore 
be available in the cross-power spectra, particularly about small large-scale signals that are correlated between redshift 
bins. The effect of different frequency window functions is shown in Fig. [5) Here the baryon oscillations only show 
up when the window is wide enough to damp down the large smaller-scale fluctuations so that the power on baryon 
oscillation scales is not dominated by contributions from smaller scales. When narrow frequency windows are used 



2 Note that at the percent level it is important to use the baryon rather than CDM velocity when calculating the redshift distortions. 

3 Note that although the CMB dipole signal has only a small effect on the dark-age 21cm power spectrum, it may make a larger contribution 
to the correlation with other sources, for example the cross-correlation with the CMB temperature during reionization (c.f. Ref. [30l0 . 
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FIG. 6: Auto power spectra for the various terms in the large-scale 21cm power spectrum. Left: narrow window function at 
z = 50 with Av = 0.01MHz (T b = -26mK); Right: broad window function at z = 40 with Av = 5MHz (T b = -15mK, o z ~ 6). 
Terms are calculated in the synchronous gauge, and 'extra' includes all effects not included in other curves; standard redshift 
distortions are defined here by the second term in Eq. (|38fl . The late-time curve is the ISW contribution from line-of-sight 
redshifting. The reionization curve is the result from sources at reionization, the other curves include the main e _2Tc damping 
effect. The standard calculation includes only monopole and redshift distortion terms; the difference from the full result is 
~ 1% at I < 50 growing to a few percent at low 1 for the narrow window function. For the broad window function that averages 
down small-scale power extra terms change the total by > 1% at I < 100 (growing to ~ 40% at low I). Note that the total 
spectrum is not just the sum of the other autocorrelation terms since it also includes all cross-terms. 



this information is hidden in the cross-correlation structure of the different source planes. 

The white-noise signal on large scales can be reduced by averaging over many redshift slices. Figure [5] shows the 
relative importance of the various terms on large scales when a broad redshift window function is used. Redshift 
distortions from scales smaller than the bin width are suppressed, and the large-scale white-noise monopole signal 
is reduced because of the line-of-sight averaging of small scale power. The relative importance of the additional 
terms is therefore larger. This raises the question of whether the large-scale 21cm signal can be useful, for example 
to learn about large scale power or as a source for the integrated Sachs- Wolfe effect 4 (ISW). There are two main 
problems. Firstly, a very broad window function is required to make the extra terms comparable to the monopole 
source, and residual monopole and redshift distortion signals will generally dominate. Secondly, since the dark age 
redshift shells are > 3/4 of the comoving distance to the last scattering surface, large-angle correlations, such as those 
due to Sachs- Wolfe potential redshifting, will be strongly correlated between redshift slices and correlated with the 
large-scale CMB. As an averaged source plane for the ISW the 21cm signal therefore has at least as much large-scale 
'noise' from other sources as the CMB and hence offers little extra information. 

Decorrelation with source plane separation is a powerful way to try and separate intrinsic and foreground signals 
due to the much smoother signal (as a function of frequency) expected from most foregrounds @, Hl| . Detection of 
small non-foreground cross-correlations is therefore particularly challenging. 



4 We thank Jeff Peterson for raising this question during a talk. 
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FIG. 7: The large-scale 21cm power spectrum at z = 20 and z = 30 (Av = 0.01MHz) if there were no Lyman-a sources, shock 
heating, minihaloes, or other non-linear effects. Solid lines are the full linear result, dashed lines include only monopole and 
redshift-distortion sources. The difference is dominated by the baryon-photon velocity term. 
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FIG. 8: The 21cm power spectrum at z = 50 for Av — {1, 0.1, 0.01, 0}MHz (bottom to top). Large widths suppress the 
redshift-distortion contribution and allow the baryon oscillations to be seen. All show characteristic damping due to line-of- 
sight averaging over the bin width at I > x/A%, and the effect of baryon pressure at I > 5 x 10 6 . 



VII. POLARIZATION 



A quadrupole anisotropy in the 21cm signal at reionization can generate polarization by Thomson scattering [321 ] . 
The signal is expected to be very small, but perhaps worth considering as polarization may be very useful to help 
with foreground cleaning of scalar modes. In principle the tensor mode signal can also be used as a cross-check on 
CMB temperature and polarization constraints on models of inflation. 
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FIG. 9: The 21cm polarization power spectra from z = 50 with Av = 0.01MHz. Thick lines show the adiabatic scalar-mode 
signal, thin lines from r = At/A s =0.1 scale-invariant tensor modes. For comparison, the instrinsic CMB polarization signals 
(labelled blackbody) are also plotted: the scalar EE blackbody curve is the lensed E-mode polarization from scalar modes 
and the scalar BB blackbody curve is the lens-induced B-mode power; the tensor blackbody curves are the CMB polarization 
power spectra from gravity waves. For both intrinsic CMB and 21cm fluctuations, E-mode polarization is shown with solid or 
dashed lines and B-modes with dash-dotted lines. There are no 21cm scalar B-modes within our approximations. Reionization 
is assumed to be fairly sharp with optical depth r = 0.09 (z la ~ 12). 



Following standard methods [33l . l34j it is straightforward to modify our scalar equations to include polarization 
in order to calculate the E-mode (gradient-like) spectrum generated at reionization. For our idealized analysis we 
neglect any effects due to magnetic fields, anisotropy of the hydrogen triplet state distribution and inhomogeneity of 
reionization. Typical numerical results are shown in Fig. [9] 

Gravitational waves (tensor modes) are known to be subdominant to the scalar modes, but can also source 
anisotropies by their anisotropic shearing. There are two mechanisms. Firstly, the metric shear can directly change 
the 21cm photon frequency along the line of sight, which distorts the emission shell in much the same was as the 
redshift distortions and line-of-sight effects do for the scale modes. Secondly, the CMB temperature at absorption 
will be anisotropic due to gravitational waves between the absorber and the last scattering surface: this causes an 
anisotropy in the 21cm absorption. At reionization the quadrupole component of these anisotropies source both E 
and B polarization. The signal is a small fraction of the blackbody tensor signal because the optical depth for 21cm 
emission is only ~ 0.02. Note that since there is no intrinsic 21cm polarization before reionization the lensing-induced 
21cm B modes are much smaller. 



VIII. NON-LINEAR EVOLUTION 

Although the dark age perturbations are quite small, non- linear effects can still be important. This is clear from 
Fig. 01 where the perturbation amplitude on Jeans' scales at z = 50 corresponds to density perturbations of order ten 
percent. We shall estimate the effect of non-linear evolution u sing Newtonian perturbation theory in the approximation 
that vorticity and decaying modes are unimportant [35L [3rl [371 [3^. [39j. The result should be quite accurate for the 
dark matter density as perturbation theory gives good results in the mildly non-linear regime. Assuming an initially 
Gaussian field one might expect the non- linear contribution to the power spectrum to be of order |A C | , corresponding 
to a correction of about one percent. However, on very small scales, there are many larger scale fc-modes that contribute 
to the local density, giving a total effect from all mode couplings that is significantly larger than the simple estimate. 

For our approximate analysis we focus on the CDM perturbations during matter domination, neglecting any effect 
due to the baryons. The first non-linear contribution to the power spectrum comes from two terms, V22 coming from 
the square of the second-order perturbation, and V13 coming from the cross-term between the first and third order 
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FIG. 11: The linear 21cm angular power spectrum at z — 30 and z = 50 (solid lines), and the approximate higher-order 
correction from non-linear evolution (dashed lines). Window functions are taken to be sharp (Au — 0). 
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FIG. 12: The 21cm power spectrum for bins at z = 50 and z = 52 (Az/ = 0.1MHz), with lensing (thin lines) and without 
lensing (thick lines). The left figure shows Ci(z,z), the right figure shows the small but non-negligible cross-correlation. 



perturbations. We give the results for the CDM density, velocity and cross-correlation power spectra in Appendix [Fj 
See Fig [10] for typical numerical results. 

To calculate the effect of non- linear evolution on the 21cm power spectrum one should really evolve the full coupled 
baryon, temperature and CDM equations to third order. This is beyond the scope of this paper. Here we instead 
make the approximation that the 21cm monopole source and baryon velocity power spectra have the same fractional 
contribution from non-linear evolution as do the CDM power spectra. Using Eq. (|39|) we then get the results shown 
in Fig.[TTJ Nonlinear effects are important at the few-percent level on small scales even at redshift z ~ 50. At redshift 
z ~ 30 there is a ~ 10% correction to the power spectrum at I ~ 10 5 . A more accurate analysis at redshift z ~ 30 may 
need to account for gas shock heating, minihaloes, or even rare first sources (see e.g. Ref. [3| and references therein). 



IX. LENSING 



The line emission spectrum is not affected by lensing at first order. However perturbations along the line of sight 
can produce a non-negligible higher- order effect. The effect on the two- and three-dimensional power spectra has 
been analysed in detail in Ref. [4fJ. Lensing acts rather like a convolution of the angular power spectrum with the 
deflection-angle power spectrum. Since the 21cm power spectrum is rather smooth, the effect is significantly smaller 
than on the CMB. Nonetheless the baryon wiggles can be smoothed at well above the cosmic variance level. Here we 
note a simple, non-perturbative approach to calculating the effect on the angular power spectrum. This is based on 
the lensed correlation function method often applied to the CMB temperature and polarization (4L|, EH • 

Since the lensing field is nearly linear the effect is well described in terms of the lensing deflection angle power 
spectrum Cf(z,z'). In linear theory this is given in terms of the primordial power spectrum 'Pfi(k) and transfer 
functions T 9 (k; rj) for the Weyl potential, W = (0 + ip)/2, by (see e.g. Ref. [43] ) 



C?(z,z') = 16ttZ(Z + 1)x 



X(2) 



X(z) - X 
X(z)x 



xO') 



dx'Ta,(k;r)o - x')ji(kx) 



X(s') - X' 
X{z')x' 



(44) 



where x( z ) is the conformal distance to the window at redshift z. We assume the window function is narrow compared 
to the distance so that the lensed sources can be approximated as a single source plane. For nearby bins the terms in 
square brackets are approximately equal. 

In terms of C"(z, z') and the unlensed spectrum Ci(z, z'), the lensed power spectrum is given to very good accuracy 
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by 

C,(z,z')«E^T^(^') £dco S pd l m (P)e- l(l+1 ^ )/2 J2 In[l(l + l)C gl , 2 (P)/2]d l n _ n ([3), (45) 



n=— I 

where I n is a modified Bessel function, d l mn are reduced Wigner functions, cr 2 (/3) = C g i(0) — C g i(/3), and 

9/ I 1 

C gl , 2 (/3) S ^^Cr(z,z')rf-ii(/3), 
i 

s E -ir C "^ (46) 
2 

The sum over n is dominated by the lowest terms with \n\ « 0. For a Au — 1MHz window at z — 50, the lensing- 
induccd smoothing of the baryon oscillations can be nearly a percent, but elsewhere the effect on the angular power 
spectrum is generally very small because it is smooth on the scale of the width of the deflection angle power spectrum 
(Al ~ 100). Figure [T2l shows the significant smoothing effect on the cross-correlation power spectrum on the scale of 
the baryon oscillations. 

In addition to the small effect on the power spectrum calculated here and in Ref. po| . lensing also makes the 
distribution non-Gaussian. Combined with tomographic information this can be used for lensing reconstruction [44], 

51, SI. 



X. OTHER NON-LINEAR EFFECTS 



On small scales much larger than the line width the non-linear angular distribution can be obtained from integrating 
Eq. HU). In the Newtonian approximation, neglecting self-absorption and lensing, and dropping small n • v terms 
that are unimportant on small scales, we have the approximate non-linear result 



/(n, z) 



(1 + A s )p s 



+ (n • dw/d x )ln 



(47) 



x' 



where \' are displaced positions satisfying \' — x( z ) ~ n ' v (x')/^- in the linear approximation we used A s (x') = 
A s (x) , however this will not be a good approximation if the Doppler displacement is comparable to the scale of the 
perturbation. The RMS velocity is around 3 x 10~ 4 at redshift 50, corresponding to a displacement of ~ 0.3Mpc, 
and hence suggesting that non-linear Doppler displacement effects may become large at k ~ lOMpc -1 . However the 
small-scale power spectrum will be unchanged under small bulk radial displacements: the spectrum at wavenumber k is 
only really sensitive to the difference in the radial displacement over a distance comparable to 1/fc. Since the velocity 
power spectrum falls with scale the effect is much smaller than indicated by the above estimate. The displacement 
on scale k is given by Sx ~ n • v/H ~ A c /k, and thus k$x will only become large on scales where A c does. Coupling 
to smaller scales will give an effective velocity dispersion similar to a scale-dependent thermal line broadening, and 
will tend to reduce the power. Since the velocity power falls with k this is never a large effect. The total non-linear 
effect is therefore likely to be comparable to the estimated contribution from just the non-linear evolution. 

A perturbative analysis using third-order Newtonian perturbation theory is given in Ref. (47| . but a detailed 
investigation of the effect on the small scale 21cm power spectrum is beyond the scope of this paper. Since higher 
order corrections to A s (x') are large on small scales, but the correction to the power spectrum is small due to the large 
coherence length of the displacement, a perturbative treatment will involve delicate cancellations between different 
independently large terms. 



XI. CONCLUSIONS 



The 21cm signal from the dark ages is potentially a powerful probe of cosmology. We have derived full linear-theory 
results for the angular power spectra, and confirmed that, for most purposes, standard approximations including 
only monopole sources and redshift distortions are accurate at the percent level. However additional velocity and 
potential terms can be non-negligible on very large scales, and optical depth and ionization fraction perturbations 
have a significant percent-level effect on all scales. Non-linear evolution can be a many-percent effect on small scales 
even at high redshift, and gravitational lensing can also be important, though this is easily modelled. 
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In this paper we have focussed on theoretical issues. The observational problems are formidable, requiring many 
square-kilometers of collection area, low radio-noise, and, for high redshifts, getting outside the atmosphere to avoid 
the opaque ionosphere; see Refs. 0, H, 0, IH[ for further discussion. The additional hurdles to measure the distinctive 
polarization signal that we calculated, which is below the intrinsic blackbody signal, are quite possibly insurmountable. 

The 21cm spectrum is much more sensitive to the residual ionization fraction after recombination than the CMB 
temperature: the ionization fraction governs the evolution of the matter temperature which in turn affects the 
background brightness and the evolution of the perturbations. Current modelling of recombination is probably 
insufficiently detailed to be able to compute the 21cm power spectrum to an accuracy of more than a couple of percent 
as there are a wealth of subtle effects that need to be modelled to calculate hydrogen recombination accurately [49| ■ As 
an example of current uncertainties, using the result of Ref. [5(j rather than RECFAST [22j for the post-recombination 
ionization fraction changes the 21cm brightness temperature at the two-percent level. If future observations are to 
achieve the ~ 0(1CT 5 ) accuracy on cosmological parameters that is possible in principle, a very detailed analysis will 
be required of the full electron, atomic state and velocity distributions, with all interactions, extending the work of 
Ref. [ljj]. Furthermore, even at redshift of 50 the Jeans'-scale density perturbations are already O(0. 1) and non-linear 
corrections are important at the many percent level. We estimated this effect with an approximate second-order 
CDM perturbation theory calculation, however accurate results will require a detailed study of the full baryon-CDM 
evolution at second and third order or beyond. 

Although we have concentrated on the power spectrum from the dark ages, much of our formalism can be adapted 
straightforwardly to the reionization epoch given models for the background evolution and additional sources. Our 
numerical code is publicly available 5 . 
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APPENDIX A: BARYON PERTURBATION EVOLUTION 

After recombination, Compton scattering has no significant effect on the baryon density evolution and the fractional 
synchronous gauge baryon perturbation Af, obeys the equation 

« \n 2 A m , (Ai) 

where (pf, -I- p c ) A m = p c A c + ptAi, and the second line assumes matter domination. We make the approximation that 
Pc 3> Pb so the baryons have no effect on the CDM evolution and we can use the usual result in matter domination 
that A m (rj) = A c (r/) = A c (ri^)ri 2 / r/ 2 . On small scales before recombination the baryons are tightly coupled to the 
photons and silk damping erases perturbations so that at recombination (77 ~ 77*) we have A;, (77*) ~ and At (77*) ~ 0. 
On scales where baryon pressure is irrelevant, kc s /H. ^ 1 at all times, the solution after recombination is 

A b /A c = l- 3 (^) 2 + 2(|) 3 . (A2) 

On smaller scale the baryon pressure is important via the k 2 c 2 term in Eq. (|A1[) . The sound speed decays approx- 
imately as I/77 while Compton scattering couples the gas temperature to the CMB temperature. Once coupling 
becomes ineffective at 77 ~ r\ a the gas cools adiabatically and c s decays as l/?7 2 . We can solve analytically for the 
evolution of the baryon perturbation for both limiting behaviours of the sound speed. For 77* < 77 < ?7 a the result is 



At/Ac = -=— 



( cosIV^ 2 - l/41n(77 / 77,)] + 5 sin[Vx 2 - 1/4111(77/77,)] 
\VJ \ 2y/x 2 - 1/4 



(A3) 



5 http://camb.info/souces 



22 



where we defined x = kc s (r]^)ri t and assumed x > 1/2. The analytic continuation of this result to x = agrees with 
Eq. (|A2[) . A similar result is given in Ref. [51J , which also discusses the generalization to non- negligible baryon fraction. 
On large scales where modes are outside the baryon sound horizon at recombination (x <C 1) the perturbations fall 
into the CDM potential wells in the order of a Hubble time and Af, — * A c . On small scales the baryon pressure causes 
the perturbations to oscillate once they reach pressure support. The relative amplitude of the oscillations about the 
midpoint decays as (?7*/f?) 5 / 2 , so by the time of adiabatic cooling we may neglect this term to ~ 10% accuracy. The 
adiabatic cooling equation is a forced harmonic oscillator equation in X/f] with solution involving sine and cosine 
integrals that can be written (for rj > rj a ) as [52j 



Aft = A cos(w) + B sin(u) + A c 



dt 



te 



i 2 + l 



(A4) 



where u = xrj a /rj — kc s (r])r) and constants A and B are defined by the initial conditions. The term in square brackets 
monotonically decreases from one as k (and hence it) increases, and describes the main effect of the pressure. Smoothly 
matching to Eq. (|A3|) (with the oscillation dropped) assuming a sharp transition at rj a , the result can be written as 



A fc /A c 



1 



cos t 
t + u 



-dt 



6 



cos(u 



■ sin(u 



(A5) 



As a function of k this solution reproduces the qualitative fall-off in power show in Fig.dJ though the decaying nature 
of the solution is well hidden in Eq. (|A5|) . The small-scale asymptotic form for k — > oo (u, x ^> 1) is 



A b /A c = \ - 12^r sin(a; - u) 



G 



[kc s (r))r]} 2 



1-0 



1< 

x rj A 



(A6) 
(A7) 



As expected the transfer function falls off as ~ 7i 2 /(fc 2 c 2 )| e on small scales with suppressed oscillations. This power- 
law fall off will only hold on scales where baryon diffusion can be neglected; on ultra-small scales the spectrum will 
be exponentially damped. Baryon pressure becomes important at fc pre ssure ~ V( c s(?7*)??*)j diffusion damping will 

1/2 

become important at ^diffusion ~ fcpressure^/^coiOmm wnere ^coii is the atomic collision time, which is always much 
smaller than r\ (see Fig. [1]). 

At redshifts of 30 to 100 residual Compton coupling actually gives c s scaling as r]~ 19 to rj^ 1 ' 6 , and approximating the 
transition in cooling behaviour as sharp may be expected to be a poor approximation on scales such that kc s {rj a )rj a 3> 1. 
In this limit, the following approximation can be used: 



A b (77)/A c (r7) 



_6_ r> rf_ 

V 3 J v , ky/ cjrjjcjrf) 



sin 



(A8) 



This is derived from the WKB solutions of the homogeneous equation with Greens method. It is an exact solution for 
all k when the gas is adiabatically cooling, but only holds for fec s (7j*)7j* 3> 1 when Compton heating is still important. 
It can be shown that this result agrees with Eqs. (|A3|) and (|A5|) in the limit x 3> 1 (and for a sharp transition to 
adiabatic cooling). 



APPENDIX B: GENERAL GAUGE AND NUMERICAL CALCULATION 

We perform numerical calculations in the synchronous gauge. For convenience we give the equations in a general 
jauge here. The multipole equations in the 1+3 conventions of Refs. [H, [54| are 



Fi 



[{I + l)F l+1 - + [S 2 i^ka - S ol h - Su±kA\ edj 



= ap s 6(e/a - E 2 \) 



5 W (A S - A)- 



Sn 



T s - T 7 3 



t [F t + S W (A HI - A Ts - A)f] - 8 lx V - (d e [ef] + fedj) . (Bl) 
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In the Newtonian gauge the acceleration A = — -0, the scale factor perturbation h = —(f), and the shear a = 0. In the 
synchronous gauge A = 0, h = h s /6 and the shear a = (h s + Qfj s )/2k, where h s and rj s are the usual synchronous 
gauge quantities. Equation (jBljl is valid in a general gauge but the individual terms are not gauge- invariant. In 
particular, under a change of velocity field, u a i— > u a + w a , the multipoles transform as 



Fn i-> F n - 



Fi h-> Fl 
F ^ f 

For completeness, the other variables in Eq. (|B1[) transform as 



I + Tied J 

wed J 



3 

(I>1) 



(B2) 

(B3) 
(B4) 



.4 



h — —(ufH + whi) + —kw, 
k 3 



A- 



w, 
w 

' ~k 



Hw 
~k~' 



Pa k ' 



(B5) 
(B6) 
(B7) 

(B8) 
(B9) 



With these relations, one can establish the gauge-independence of Eq. (|B1 
The line-of-sight solution is 



F = e~ T < I d V e-^ +f \ \{ap s (A s -A)- af s (A HI - A Ts - A) J) S(e/a - E 2X ) + hedj - f c F ] ji(k X ) 

T 

{Ps ~ r s f)d e [ed(e/a - E 2 i)} v - p S l ^ ( v ~ v 7 )6(e/a- E 2 i) 



+ (kA - t c v)edJjl(k X ) - ap s S(e/a - E 21 ) 7fr ^ r £(2i' + l)Q v i 1 ' ft, ( --— ) 



7 /' = 



Z'=2 



(BIO) 



(for I > 1), the spin temperature perturbation is 



A Ts - A T7 = (Ry - Rg)SC w + C w {R g A Tg - F 7 A Tt ) + -r^ 10 Ci 

Z 1 tIt 



Tg — T 



A HI -^-^-A + A Cl0 + 2A Tt (C 10 F 7 - 1) - 2R^6C 10 + ^ ¥ 



An 



A T J 



(Bll) 



and the gas temperature perturbations evolve with 



a 2 1, oU &aa T p~ ( x e 

At = — kv — 2h ■ — - 

3 3m e c(l + /ho + x e ) 



A 



l + x e /(l + / He ) 



^(A Ts -A T7 ) 



(B12) 



where we again neglected helium fraction perturbations. Equation (|35[) in the main text holds in any gauge. 

We choose to work in the synchronous gauge, evaluating the line-of-sight solution integrated over a frequency 
window function. We use the synchronous gauge because this is stable for isocurvature mode evolution, and because 
conventional calculations invariably use the baryon (or CDM) power spectrum in the synchronous gauge as output by 
CMBFAST or CAMB. Since both v — v c and t v are small on scales where our calculation is applicable (much larger 
than the line width), in the synchronous gauge [v c = 0) terms involving t^v can be dropped to good accuracy. 

We define a window function Wf{e) so that we observe J deWf(e)f. Then integrating over energies, assuming the 
window function is much broader than the line width, we use the functions 



dee 7 



'W f {e)ap s d{e/a-E 2 



1 - e~ r " 



a 2 p s W f (e v ) 



(B13) 
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G T {ri) = / dee f - T 'W f (e)afJS(e/a-E 21 ) = 



1 _ T ne " 



1-e- 



(B14) 



7fa) = / de e ^W / (e)e9e/ = -e- T "^W / ( e ,) + J tof*tf)-^ 



9 /l- e V a{r)')p a {r)') 



H{rf) 



W f (e n >), (B15) 



where W/(e) = — Wf(z)E2i/e 2 = W/(a)/£ , 2i- For 21cm we assume we measure some averaged brightness temperature 



n 



dvW T (v)T b (v) 



(B16) 



Then 



ch 3 

f h = I daW T (a)^ef 



where Wt{o) is a window over dimensionless frequency, so 



W f (e) = ^oW T (a) 



(B17) 



(B18) 



APPENDIX C: EVALUATION OF ACo 

The contribution from radiative transitions in the 21cm line to the evolution of the spin temperature is 



Nr. 

(l+N v0 )Aw-3-^A w Ko 

iVi 



4p.A w ( l-jAf v o 



(CI) 



Here, recall, ACo is the isotropic part (in the gas rest frame) of the photon occupation number integrated over the 
line profile. The occupation number has contributions from the background and perturbed CMB and the 21-cm line 
radiation, i.e. 



K= kBT 1{1 + Q) + ^ if+Sfy 

making the Rayleigh- Jeans approximation for the CMB. Evaluating 

1 



(C2) 



A4o 



An 



dEdVtM{E : - E 2 i) 



1 f E 
= — / dEdQ ^Af(E,e)<f>(E - #21 ) 

47T J E 



(C3) 



for the CMB, we find 



A'' 



(CMB) 



1/0 



i 



(i 



(C4) 



The 21cm part takes more work. First consider the contribution from 5f; since Sf is first-order, we only require 
its monopole in the conformal Newtonian gauge. As we are integrating over the line profile, at any time r\ we are 
including only 21-cm radiation that was 'produced' at earlier times within a narrow time window A77 ~ AE / {TLE21) 
where AE is the line width. We therefore generalize Eq. (fT8|) for e very close to e v under the assumption that the 
perturbations have (comoving) wavelength ^> Arj. Extracting the monopole, we find 



-Z5fo(r},x, e) 



_L U _ I V • v ) f e- f + (A s + V) (1 - e- f ) - (A Hl - A Ts + [l - (1 + f ) e - f ] 



(C5) 
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The divergence of the gas velocity arises here from the monopole of the redshift-space distortion: 

1 f 1 
— / dQ e e- (e- Vv) = - V • v. 

Integrating over the line profile, we obtain the following contribution to A/^o from 21cm perturbations: 



(C6) 
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(A s + i>) l 



1 - e~ r " 



(Am - A Ta + VO 1 + e- T " - 2 



1 - e- r " 



1 



-V-v 



1 - e- T » 



We also need the contribution from /. Correct to linear-order in v, this is simply 



(/) 



v 



1 - e- T " 



(C7) 



(C8) 



non-linear corrections in v may be important if i?2i|v|>A_E. Combining all results, and expanding to first-order in 



t v , we obtain our desired result 



l + A s + ^ + l U-iv-v 



(C9) 



In this approximation, the perturbed Rayleigh-Jeans' brightness temperature, T*N v q, that we use in the main text is 



T 7 + T b0 = f 7 (l + A T7 ) + -1(T, - T 7 



1 + A s + ip 



H 



-V • v 



(CIO) 



As we might have anticipated, this can be expressed in terms of the (perturbed) optical depth t v of Eq. [30l and the 
monopole of the CMB temperature T 7 as T^o = t v (T s — T 7 )/2. 



APPENDIX D: IONIZATION FRACTION PERTURBATIONS 



For our approximate analysis of ionization fraction perturbations we start the perturbation evolution after Helium 
has recombined, so we take rim = (1 ~ x e) n H an d n e = x e nn where Uh is the total number density of ionized and 
unionized hydrogen. We then use the effective equation of RECFAST [22| 



x e = -aC r (x 2 e n H a - /3(1 - x e )e~ EH2 ° /kBT <>^ 

where Eh2s is the energy of the transition from the ground to the 2s state and 

= 1 + JCAn H (l - x e ) 
r = \+K{K + P)n H {l-x e )' 

The recombination and photoionization coefficients are related by 



' m e k B T g 



3/2 



-E 2s /k B T g 



where E^ s is the ionization energy from the 2s state, and a is also a function of the temperature fit by 

a = F- 



a a (T g /\Q A K) b 3 _ x 



l + c{T g /lO A K) d 



m s 



(Dl) 



(D2) 



(D3) 



(D4) 



Here F is a fudge factor taken to be 1.14 and a a = 4.309 x 10~ 19 , b = -0.6166, c = 0.6703, d = 0.5300. The constant 
two photon 2s-ls decay rate is A = 8.22458s -1 . Dependence on the expansion rate enters through the cosmological 
redshifting term /C, given in the background by JC = a\^ I2p /(8nTi.). 

A perturbed version of Eq. (|Dip may be inappropriate due to the effect of perturbation velocities on the escape 
probabilities that went in to deriving the result. However the detailed perturbation evolution during recombination 
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does not affect the 2fcm absorption signal as long as the residual perturbations after recombination are correct, at 
which point recombination is limited solely by the low rate of electron capture to an excited state. A full analysis is 
beyond the scope of this paper, so we proceed from the perturbed version of Eq. (|Dlj) and argue this is sufficient: 

A Xe = — (A Cr +A a - A Xe - A) 

- aC r | (2A Xc + A H )ax e n H + j A Xs - Q + ( ±- - lj A Tg } p e - E *"/ k » T ' j . (D5) 

We generalize JC to a perturbed universe as IC — 3X H2p / (8nV a Ug) to account for the local baryon expansion rate; this 
should be correct on scales sufficiently large compared to the mean Lyman-a photon interaction length (which is very 
short). The perturbed terms are then 

. plCrm {[A H + A K + A^(l + EAn g (l - x e ))} (1 - x e ) - x e A x J 

Cr ' (l+lCAn H (l-x e ))(l+lC(f3 + A)n„(l-x e )) [ ' 



3 E- 



2* 



2 k B T t 



(D7) 



b + c(T g /lO^KY(b-d) 

Aa = l + c(T g /WKY At « (D8) 

h kv , . 

= ~H~ SH~ (D9) 

In the synchronous gauge the dominant term on large scales is initially due to fluctuations in the expansion rate, 
Ax; ~ At,/37i: A Xr is the same sign as A c because overdensities expand less fast and hence have a higher ionization 
fraction than the average (A^ appears to have been ignored in e.g. Refs. [H^,[5^]). At later times when i e «l the 
main effect is from perturbations in the hydrogen density and temperature; this leads to A Xc being the opposite sign 
to A c at late times because overdensities recombine more efficiently. Note that at late times some terms in Eq. (|D5[) 
are invalid because At Ei/ksTg becomes of order unity, however the error is harmless because the entire term is 
exponentially suppressed: at late times the evolution is given approximately by 

x e « —aux\nii (D10) 
A Xe sa -aax e n H (A Q + A Xe + A H ) ■ (Dll) 

It is not necessary to model the early evolution correctly to get approximately the correct late-time answer. Indeed the 
late-time evolution neglecting velocity effects should be quite accurate as recombination is limited by the low electron 
capture probability. However the overall ionization fraction evolution is limited by the precision of the RECFAST 
model, in which a single fudge-factor accounts for deviations of an effective three level atom model from the full result. 

Using the equations given here the effect on the 21cm power spectrum of neglecting A Xis is 0(2%) at z ~ 50, 
almost entirely due to the indirect effect on the evolution of the temperature perturbation. Note that, at the Jeans' 
scale, ionization fraction perturbations also have a linear effect on the baryon density evolution due to the modified 
evolution of the gas temperature perturbation (and hence the baryon pressure perturbation). 

APPENDIX E: SPHERICAL BESSEL FUNCTION INTEGRALS AND APPROXIMATIONS 

Integrals of products of spherical Bessel functions can be done analytically using the result 

Trr (n + i)r( '+' , + 1 -" ) 

drji(r)ji'(r)r = ^ n+2 ^ ^_;/ +2+7t ^ p ^ W / +3+7 ^ p ^/_; +2+ ^ y- ( E1 ) 

In particular, using this result in combination with the spherical Bessel equation we have 

f°° 1 
1(1 + 1) / dlnr^(r)] 2 - - (E2) 

Kl + 1) j d In rji{r)j\'{r) = -- (E3) 
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where the approximation holds for high I of most interest in this paper. 

At high I the Bessel functions oscillate very rapidly compared to the scale of variation in the power spectra. Direct 
numerical integration becomes slow, but quite accurate results can be obtained by averaging over the oscillations. 
Defining v = I + 1/2, for r 2 jv 2 > 1 + ^~ 2 / 3 we can use the approximation 57 | 
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Neglecting oscillatory parts that closely average to zero we then have for I 3> 1 
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APPENDIX F: NONLINEAR CDM POWER SPECTRA 



At a given redshift the two leading corrections to the power spectrum are given in terms of the linear-theory matter 
power spectrum at that redshift, V{k) = k 3 P(k)/(2ir 2 ), by [H, [H, [H, i| 
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For the small scales of interest to us here the dominant contribution comes from q < k, with significant mode coupling 
to all scales where the spectrum is growing logarithmically. This can become numerically difficult because for r<C 1 
the two terms P22 and Pi 3 become large but almost cancel. For r < r s we therefore use an approximate series 
expansion, switching to the full result at r > r s . The expression for the correction to the matter power spectrum can 
then be written 
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where r s < e < 1/2 (numerically it is better to chose e ~ 1/2). For a scale invariant primordial spectrum the small 
scale power spectrum is V{kr) ~ P(/c)[l — ln(r)/ ln(ro)] 2 , where rg corresponds to the much larger scale fco where 
logarithmic growth starts. The main non-linear contribution to the power spectrum can then be approximated from 
the first term in Eq. (|F3j) as ~ — ln(ro)[P(fc)] 2 ~ \n(k/ka)[P(k)] 2 . On the Jeans' scale at z = 50 this implies a 
fractional second-order contribution to the power spectrum of ~ O(10)V(k): on small scales non-linear effects are 
more important than one might naively think [37| . See Fig [10] for typical numerical results. 

The 21cm angular power spectrum on small scales is also sensitive to redshift distortions, for which we need to 
estimate the second-order velocity power spectrum and the cross-correlation with the density. The results are similar 
to those for the density with 
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where "P™ is the power spectrum of kv/Tt. The series result for use at high I is 
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Similarly the cross-correlation power spectrum is given by 
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with the series result 
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